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In random heteropolymer globules with short-range interactions between the monomers, 
freezing takes place at the microscopic length scale only, and can be described by a 1- 
step replica symmetry breaking. The fact that the long-range Coulomb interaction has no 
intrinsic length scale suggests that freezing in random polyampholyte globules might take 
place at all length scales, corresponding to an overlap parameter q(x) that increases 
continuously from zero to its maximum value. Study of the polyamphol5^e globule within 
the independent interaction approximation seems to confirm this scenario. However, the 
independent interaction model has an important deficiency: it cannot account for self- 
screening, and we show that the model is only reliable at length scales shorter than the 
self- screening length. Using the more realistic sequence model we prove that in the 
general case of a random heteropolymer globule containing two types of monomers such 
that unlike monomers attract each other, freezing at arbitrarily large length scales is not 
possible. For polyampholyte globules this implies that beyond the self- screening length, 
the freezing behavior is qualitatively the same as in the case of short-range interactions. 
We find that if the polyamphol5^e globule is not maximally compact, the degree of 
frustration is insufficient to obtain freezing. 

I. INTRODUCTION 

A polyampholyte is a polymer containing both positively and negatively charged 
monomers. Despite the fact that polyampholytes provide simple models for electrostatic 
interactions in biopolymers such as proteins, they have been studied intensively only 
during the last 10 years. Most attention has been paid to the conformations of isolated 
molecules in solution. ^"^ The authors of Ref. 1 considered a chain with zero net charge, 
containing a fraction / of charged monomers randomly distributed along the chain. 
They showed that if the chain is sufficiently long, it collapses to form a globule (see also 
Ref. 7, where it was shown that oppositely charged polyelectrolyte chains merge to form 
globular complexes). Generally, the globular state^'^ of a polymer chain is characterized 
by a well-defined, non-zero density, and by the fact that its correlations are short-range 
(short compared to the linear size of the globule). Thus, the globular state differs 
qualitatively from the coil- state. Depending on the value of / , this globule may or may 
not be maximally compact. In either case, further collapse is prevented by excluded 
volume interactions.^ In order to find the density of non-compact polyamphol5l;e 
globules, the authors of Ref. 1 used the Debye-Hiickel theory^^ to describe the charge 
density fluctuations inside the globule. This is justified provided that the fraction / of 
charged monomers is small, and the quenched sequence of charges is random.^ Things 
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might change if the polymer carries a net charge. In Refs. 3,4 it was shown via computer 
simulations that if each monomer carries a charge ± e , then even a small net charge of 
order e-Z/V is sufficient to stretch short chains (A'^ is the number of monomers per 
chain). In Ref. 5 the effect of a net charge was investigated theoretically. It was shown 
that there exists a critical value oc ^f]^ for the net charge such that polymers with 
1 2n ' < 2* collapse into spherical globules, whereas polymers with 1 2n ' > 2* collapse 
into elongated globules. The temperature dependence of the characteristics of the globule 
was derived in Ref. 6 using a Flory-type theory. Multi-chain effects were studied in Ref. 
11. It was shown that for neutral polyampholyte samples in solution the globules in the 
supernatant are spherical, since molecules with an opposite excess charge form neutral 
complexes (the same as what happens in mixtures of oppositely charged 
polyelectrolytes^) . 

The aim of this paper is to investigate the possibility of "freezing" in polyampholyte 
globules, and in particular we are interested in the influence of the long-range character 
of Coulomb interactions. The freezing-phenomenon has already been studied extensively 
for the situation where the interactions between the monomers are short-range. ^^"^^ First 
we give a brief explanation of the concept of freezing in compact heteropolymer systems. 
Consider a heteropolymer consisting of several different monomer species, which are 
distributed randomly along the chain. At low temperatures the (compact) system tries to 
reach a low-energy state. However, due to the constraints imposed by the chain 
connectivity, it is not possible to minimize all interaction energies at the same time. 
Therefore, the system is "frustrated," and has many local free energy minima. At low 
temperatures the ergodicity is broken: the phase space of the polymer globule is divided 
in many phases (physical states). As in spin glasses, these phases are not related to each 
other by symmetries of the Hamiltonian (see Ref. 19 for a review on spin glasses and 
replica theory). It was conjectured in Ref. 12 that to a first approximation the energy 
levels of a polymer globule are independent from each other, and that they have a 
Gaussian distribution. This would imply that a compact heteropolymer globule with 
short-range interactions is equivalent to the random energy model, as far as the 
thermodynamics is concerned. This was shown to be true for a microscopic model in Ref. 
14. The ergodicity breaking of the random energy model can be described by a 1-step 
replica symmetry breaking. The physical meaning of a 1-step replica symmetry 
breaking for polymers is as follows. Consider n » 1 heteropolymer chains, all having 
exactly the same monomer sequence (such chains are called replicas). At conditions 
where ergodicity breaking occurs, the conformations of two replicas either coincide 
exactly up to the microscopic length scale (modulo rotations and translations), or they 
have no structural relationships at all.^"* It follows that when the ergodicity is broken, the 
relative monomer positions within one chain are more or less fixed. One describes the 
situation by saying that the chain is "frozen" in one conformation. Protein folding is a 
special example of heteropolymer freezing, where the energy of one of the frozen 
conformations, called the native one, is much lower than that of the others.^^ The 
applicability of the random energy model to describe ergodicity breaking in polymers is 
not unlimited. For instance, it was shown in Ref. 15 that for low space dimensions d <2 
the replica symmetry breaking is of a more complicated nature. The reason for this 
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qualitative difference between globules in J > 2 and d <2 space dimensions can be 
explained in a simple way. As an example for the case d >2 consider a globule in 3 
dimensions. Let r be a point visited by the chain, and consider the sphere with center r 

and radius R . The random walk describing the chain conformation will need roughly R 
steps in order to leave this sphere, giving it an average density p oc XjR . Therefore, a lot 
of empty space is left, and the chain has to return many times in order to create a compact 
globular state. The chain returns because it is "reflected" at the boundary of the globule, 

and relatively short pieces of chain with length oc n"^'^ traverse the whole volume. It 
follows that the chain is rather self-entangled, and so it is difficult to change the structure 
of part of the globule without affecting the structure of the remaining part. Therefore, two 
conformations that are different in some region of space will probably be different 
throughout the globule. This structural independence of different conformations makes 
that the random energy model provides a good description for d = 3. However, this 
reasoning breaks down for d = 2 , because then the number of random walk steps needed 
to leave a circle with radius R is proportional to the surface of the circle. In other words, 
if the polymer chain reaches the point r , it will completely fill up the neighborhood of 
this point before moving on to other regions. This leads to a so-called crumpled globule, 
in which monomers that are close to each other in space are also close to each other along 
the chain. It follows that globules m d <2 dimensions are always crumpled. Since the 
various parts of a crumpled globule do not interpenetrate and are therefore independent 
from each other, it is conceivable that the conformations of two crumpled globules 
coincide in some region of space, while they are completely different outside this region. 
Consequently, the assumption of structural independence of different conformations, the 
cornerstone of the random energy model, breaks down in low dimensional space. Another 
vital ingredient for the applicability of the random energy model is the short-range 
character of the interactions. If the interaction has a length scale L larger than the 
microscopic length scale, then by performing a renormalization transformation one can 
see that the random energy model is applicable on length scales R^L , but it is not clear 
what happens at shorter length scales. If the interaction has no characteristic length scale 

at all (as it is the case for the Coulomb potential V(r) oc r^~'^ ), the random energy model 
might not be applicable at any length scale. Application of the qualitative investigation of 
Section 4 in Ref. 15 to polyampholytes seems to confirm this, and suggests the existence 
of a critical dimension d = 6 below which the random energy model breaks down, 
implying a non-trivial replica symmetry breaking. However, as we will see, part of this 
conclusion is an artefact of the model employed (the independent interaction model). As 
suggested in Ref. 24 and confirmed by the present paper, at length scales larger than the 
polymer self-screening length the random energy assumption becomes valid, since the 
Coulomb interaction is effectively short-range at those length scales. However, it is 
worthwhile to investigate what is happening on smaller length scales. This might also be 
relevant for the problem of protein folding. Of the 20 natural occurring amino acids, two 
of them are positively charged, and two of them are negatively charged. At first sight it 
seems that Coulomb interactions are irrelevant because proteins fold under physiological 
conditions where the salt concentration is high and electrostatic interactions are 
effectively screened. However, screening can only take place in the coil-state of the 
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chain. Once the protein is collapsed into a globule, all water and salt is squeezed out, and 
the electrostatic interactions between charged groups in the interior are unscreened. 

The main part of this paper deals with freezing in the globular state of a weakly charged 
"statistically neutral" polyampholyte chain in a 9 -solvent or a poor solvent, because for 

these situations a rigorous analysis is possible. As described above, the globular state can 
arise in various situations. It can arise by the collapse of a single chain, leading to an 
elongated or spherical globule, depending on the net charge.^ The globular state can also 
arise by the simultaneous collapse of two chains whose net charges cancel each other, 
leading to a spherical globule. Finally, the globular state is present in the precipitate of a 
polyampholyte solution. Although in the latter two cases the system consists of several 
chains, the single-chain analysis of this paper remains valid provided that the chains are 
long enough. This is because even in the case of a single long globular chain, monomers 
that are far separated from each other along the chain (but that are possibly close to each 
other in space) do not "know" that they are attached to each other. This is due to the short 
correlation length inside the globule,^ making the long loop connecting these monomers 
invisible. 

Section II provides a summary of the model introduced in Ref. 1 to describe the globular 
state of a weakly charged polyampholyte. Section III presents the theoretical tools 
developed in Ref. 15 to analyze freezing in situations where the self-overlap of the 
replicas is small. In the case at hand this criterion is satisfied: if there is freezing, each 
replica will be frozen on the self- screening length scale, which is much larger than the 
microscopic length scale since the polyampholyte is only weakly charged. In Section IV 
we analyze the polyampholyte globule using the independent interaction model, leading 
to the prediction that in the absence of salt the overlap parameter^"* increases continuously 
from zero. In Section V we show that the independent interaction model provides a 
correct description only at length scales smaller than the self- screening length, and on 
larger length scales the solution found in section IV must be rejected. By using the more 
realistic sequence model we prove that the overlap parameter cannot attain arbitrarily 
small values, and so q(x) must be discontinuous even in the absence of salt. In Section 
VI, finally, we summarize our results, and discuss briefly the possibility of freezing in the 
globular state of highly charged polyamphol5l;es. 

II. PHYSICAL MODEL 

To describe the globular state of a weakly charged polyampholyte we follow Higgs and 
Joanny.^ They considered an overall neutral random polyamphol5l;e chain containing 
both charged and neutral monomers. Due to electrostatic interactions, the chain is 
collapsed to form a globule. Since the chain is weakly charged, the globule does not have 
maximum density, in other words, there is still solvent present inside the globule. We 
mainly consider the case where the solvent is either under -conditions, or under poor- 
solvent conditions with respect to the neutral monomers. In both cases the chains obey 
random walk statistics on all length scales^ in between the microscopic length scale and 
the linear size of the globule. This property is used explicitly in our calculations. If the 
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solvent is good, then at length scales larger than the blob size (which equals^ the polymer 
self- screening length) the chain statistics is again Gaussian, and on these length scales the 
behavior of the system is qualitatively the same as that under -conditions. In this paper 
it will tacitly be assumed everywhere that the solvent is under 9 -conditions, unless stated 
otherwise. The composition of the chain is determined by the fraction / « 1 of charged 
monomers. Whether the number of neutral monomers in between two charged monomers 
is constant or random is of no importance. Apart from the sign, all charged monomers 
carry the same charge ±e. The extrapolation of the results to multiply charged 
monomers is trivial. There is no charge correlation along the chain; in other words the 
sequence of charges is BernouUian (random). The typical root-mean-square distance a 

between two charged monomers neighboring along the chain is given by -b^ I f > 
where b is the size of a neutral monomer. It was shown in Ref. 1 that in the globular 
state, further collapse is prevented by the excluded volume repulsion (the entropy loss 
due to confinement is a surface effect and can be neglected for long enough chains 
(volume approximation)). In case of non-compact globules the excluded volume 
repulsion is quantified by the third order virial coefficient in case of a 6 -solvent or a poor 
solvent, and by the second order virial coefficient in case of a good solvent. Under 9 - 
conditions the resulting globule density p , defined as the number of charged monomers 
per unit of volume, is given by^ 

where i = /Anek^T is the Bjerrum length; l~Q.lnm in water at room temperature. 
We require that b> i , because otherwise the charged monomers would stick together. In 
deriving Eq. (1) the authors of Ref. 1 assumed that the charge density fluctuations inside 
the globule can be described by the Debye-Huckel theory, which is justified^ if the 

polymer self-screening length k"^ is much larger than the average distance p"^'^^ 
between two charged monomers. Using Eq. (2) one finds that this condition is equivalent 
to » f f' . At the same time, this implies that the globule density as given by Eq. (1) 

is much smaller than the maximum density ^^^^ = f jb^ . The globular state can be 
considered as a melt of blobs.^^ Both under -conditions, and under good solvent 
conditions the blob size coincides^ with the polymer self-screening length k"^. The 
number of charged monomers per blob, and the radius of a blob are given by^ 

,2 I 2 

A^b--^»1 ^b-Kli-— (2) 

fl^ b p 

The melt of blobs is moderately strong interacting, which can be seen in the following 
way. The characteristic interaction energy E between two neighboring blobs can be 

estimated as E /AtieR^ ~ k^T , where g^, ~ eJN^ is the characteristic charge of a 
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blob. Since the Coulomb interactions are screened at distances larger than , on length 
scales R » R^^ the globule is equivalent to a (maximally compact) chain with length 
N/N^ having random short-range interactions (the fact that this equivalence also holds 

with respect to the freezing behavior is one of the conclusions of this paper). In Ref. 16 it 
was shown that such a system has a freezing transition when the dispersive interaction 
energy between different monomers attains a critical value, which is of the order k^T 
(corresponding to Xpiory ~ 1 )• So, it may be concluded that the polyampholyte globule is 

at least not too far away from a freezing transition. This means that a considerable 
amount of frustration must be present in the system. The origin of this frustration can be 
understood in the following way. The physical state of a polyampholyte chain is 
dominated by conformations in which the charges screen each other (note that in random 
conformations the charges are unscreened). In order to achieve this so-called self- 
screening, the charges have to rearrange spatially. This rearrangement is hindered by the 
chain connectivity, leading to a frustrated situation. 



III. GAUSSIAN APPROXIMATION 

On length scales smaller than the blob size, the thermal energy exceeds the electrostatic 
energy, and so freezing could only occur on length scales R ^ k"^ . Since it follows from 

Ny^ » 1 that Kp^ is much larger than the microscopic length scale a , it is justified to use 

the procedure developed in Ref. 15 by Shakhnovich and Gutin. In this section we present 

a brief summary of their method. Within the framework of replica theory^^ the onset of 
freezing in a polymer globule is signaled by the appearance of a non-zero value for the 
overlap order parameter Q^^ir, r') , which is for a ^^^^ |3 defined by^"^ 

Qa^(r,7') = 1^5,(7-71) 6,(F-?pO^ (3) 

5^ is the Dirac delta function smeared out over the microscopic distance a , the vector 

is the position of the (charged) monomer in replica a, and the brackets denote a 

thermodynamic average. The order parameter Q^^^ appears in a natural way during the 

calculation of the free energy, but it is convenient to define a simpler overlap parameter^"* 
^ccp having a clear physical interpretation: 

?ap - ^ I ^ <2ap ^) = £ (Fj - Fp'- )j (4) 

The quantity between brackets gives for replicas a and P the fraction of monomers that 
occupy the same position in space (measured within accuracy a). The normalization 
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In makes q^^ dimensionless and ensures that = 1 if replicas a and [3 have 

perfect overlap up to the microscopic length scale a . Another important characteristic of 
the system is the density profile p^jir) , which is defined by 

PaW = ^^S,(?-Fj)^ (5) 

The density is in fact independent of a . In deriving the equation for the free energy we 
will need an expression for the entropy 5„ = /cg In f2„ , where n„ is the (weighted) 
volume in the phase space of n ideal non-interacting chains, associated with a certain 
density profile po((r) and overlap parameter Q^^(r,r'). Generalizing the method 

presented in Refs. 8,9,15 one can derive the following relations: 

kg J (p 

Qa^(r,r') = ^cB^d(x^-r) 6(xp-F') (pg(p (6) 

p„(F) = |i/i„5(i„-F) cp^cp 

where Jx^ = dxi--- Ji„ , and cp(Xoj) is some unknown function. The integral operator g 
is defined by^ 

8<?ixJ = ^dy^i?ixjY[8(^^-yJ (7) 

a 

g(r) = g(r) is the single bond distribution function; in other words, it is the probability 
density that in the ideal system, two (charged) monomers neighboring along the chain are 
separated by a distance r . Using the fact that the characteristic length scale R of the 
function (p (which is the scale at which freezing occurs) is much larger than the 

characteristic length scale a of the function g(r) (because R^K.~p^ » a), it is justified 
to make the approximation 

a, J " 

where j = 1, 2, 3 denotes a spatial index. With this approximation, the relations in Eq. (6) 
become 
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(9) 



In Section IV the free energy F will be written as a fiinctional of (p(Xp() . In order to find 
the mean- field expression for the overlap order parameter Q^^ , the free energy has to be 
maximized with respect to (p . Following Ref. 15, we consider Gaussian trial functions 
of the form 



}) = % exp(- { ^ k^px^ ■ ) 
cc,P 



(p2 ^ ;i3(l-.)/2p^3/2 



(10) 



The matriK k^^ is of the Parisi-type; the value k^ at the diagonal may be different 

from zero. The dot represents an inner product in 3 dimensional space. Because of 
translation invariance on length scales small compared to the linear size of the globule, (p 
must be translation invariant, which is equivalent to the condition 



Va 



(11) 



The quantity D appearing in Eq. (10) is the determinant of the matrbc obtained^^ by 
deleting the a* row and the a* column of the matrbc k^^ ; this determinant is 
independent of a. After substitution of the trial function Eq. (10) for (p into Eq. (9), one 
finds the following expressions for the density and the overlap parameter: 



2ap('i''i) = 



3/2 



exp 



(1-^2) 



D. 



3/2 



Pa(^) = P 



2ap(^) = pexp 



(12) 



The Parisi matrbc D^^ is the determinant of the matriK obtained from k^^ by deleting the 

a* row, the a* column, the p* row, and the p*^ column, divided by D .^^ The 
physical meaning of Eq. (12) is illustrated in Fig. 1: the conformations of replicas a and 
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P coincide if they are measured with accuracy R ~ -Jd^ . One should realize that Eq. 

(12) represents a class of trial functions, and that the real overlap parameter might be 
different. According to Eq. (12) a weak overlap between two replicas a and P is 
achieved by assuming that their conformations are confined within the same "tube" 
whose radius R » k"^ is constant throughout the globule. In reality, however, it is 

possible that in part of the globule the conformations coincide down to the self- screening 
length, while in the remaining part they are completely different. In fact, this is probably 
what happens in low-dimensional d <2 globules with short-range interactions. 
Nevertheless, even in this case the class Eq. (12) of Gaussian trial functions describes the 
situation qualitatively. Adopting this approximation, the expression for the entropy S 
becomes^^ 

2 1 X 

— = bm--^ = —\^K{x) K{x) = -xk{x)+\dyk{y) (13) 

n^on ks 4 X J 

The one-to-one correspondence between the function k{x) and the off-diagonal elements 
of the Parisi matrix k^^ is defined in the usual way.^^ The function K(x) is related to the 
function D(x) (corresponding to the matrix D^^ ) via^^ 



D(x) = 2 



1 

C(x) J 



dy 



xK(x) J y^K(y) 



(14) 



IV. INDEPENDENT INTERACTION MODEL 

In this section we study the weakly charged polyampholyte globule in the independent 
interaction approximation,^"^'^^ which is based on the simplifying assumption that the 
coupling constants {5^}, determining the interaction strength between charged 

monomers i and j , are independent from each other. This assumption is violated in the 
case of interacting electric charges, where Bg is proportional to the product of the 
charges carried by monomers i and j . For instance, if Bjj < and Bi^. < , then 
necessarily > . However, our primary interest is to find the influence of the long- 
range nature of the Coulomb interaction, and the model is capable of capturing this 
feature. In Section V we will drop the independent interaction assumption by studying the 
more realistic sequence model. The great advantage of the independent interaction model 
with respect to the sequence model is, that it is technically easier to handle, which is the 
reason why it was considered in the earlier theoretical papers concerned with freezing in 
copolymers with a short-range interaction. ^"^'^^ We will see that, within the Gaussian 
variational approach used throughout this article, the independent interaction model can 
be solved analytically. However, in case of Coulomb interactions a very important 
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drawback of the model is that it is not capable of describing screening, which can be seen 
in the following way. Consider any monomer, say monomer / . This monomer tends to be 
surrounded by a cloud of monomers k for which 5,^ < . Now suppose that monomer / 

is approached by a repelling monomer j (that is, B^- > 0). In the physical system of a 
mixture of positive and negative charges, monomer j is attracted by the cloud, thus the 
unfavorable interaction between / and j is screened. If, on the other hand, the interaction 
strengths are independent, then the fact that / repels j , whereas / attracts k , does not 
give any information about the interaction between j and k . Therefore, as experienced 
by j , the cloud is just a random collection of charges and the unfavorable interaction 
between i and j is unscreened. This suggests that the independent interaction model can 
only describe the polyampholyte globule correctly on length scales that are smaller than 
the polymer self- screening length k"^ , which is confirmed by the analysis in Section V. 

Since it will turn out that on scales R » k"^ the system is qualitatively equivalent to a 

globule with short-range interactions, the independent interaction model describes the 
most promising regime for finding non-trivial replica symmetry breaking. On the other 

hand, at length scales smaller than k"^ the thermal energy becomes of the same order of 

magnitude as the electrostatic energy, so there is freezing at a small length scale it has 
to be weak (see also the discussion at the end of Section II). 

We consider the physical system described in Section II, and model it using the 
independent interaction approximation. We will assume a Gaussian distribution for the 
interaction strengths By . At first sight, this choice seems difficult to justify, because after 

renormalization. By is the product of two variables with a Gaussian distribution (these 

variables represent the charges of the coarse grained monomers), which is itself not 
Gaussian. However, we note that expansion of the sequence model in the regime 

R « Kp^ leads to the same expression for the free energy as the independent interaction 

model with a Gaussian distribution for B^ (compare Eq. (23), which is valid for the 

independent interaction model, with Eq. (48), which is derived from the sequence model). 
In the continuous representation of the chain, where the parameter t "numbers" the 
charges, the partition function can be written as 



Z= f J?(Oexp ^— \dt{^^\ -- {dtdt'B(t,t')V(r(t),r(t')) 

J 2a J [ dt 2J 



(15) 



In writing down Eq. (15) it was tacitly assumed that the integration over chain 
conformations r(t) is restricted to conformations having a constant density p(r) = p as 
defined in Eq. (5) (see the discussion following Eq. (41) for a justification of our neglect 
of density fluctuations). The Gaussian probability distribution for B(t,t') has mean zero 
and variance unity. 
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P(B(t, f)) = -^cxp[-^B\t,t')\ 
V27U ^ 



(16) 



For the potential V we take the screened Coulomb potential^ 

exp[-Kr] 



V(r,r')=V(\7-r\) 



V(r). 



(17) 



Within the limitations of the independent interaction model, K ^ represents the screening 
length due to the presence of salt, in which case it is an independent parameter of the 

theory. We will see, however, that this physical interpretation for leads to 

qualitatively wrong results. When the independent interaction model is obtained from the 
sequence model by expansion of the free energy in powers of the potential (see Section 

V), turns out to be the self-screening length k"^ of the polymer globule, which 

means that it can not be considered as an independent parameter. This fact is closely 

related to the observation that charges interacting via independent interactions cannot 
screen each other. Nevertheless, it is instructive to analyze the independent interaction 

model for general values of . Following the same kind of derivation as presented in 

Refs. 16,17 one arrives at the following expression for the partition function [Z"] of n 

replicas averaged over all sequences: 

Y^l^^i^^dr^dr^ !2ap(^l>^2) Gap (^3 '^4) ^(^1.^3) ^(^2' ^4) 



'^nbap)=^Bln|^^?(Oexp -^Y^^dt 



drjt) 
dt 



5[!2ap(^>^')-eap(^,^')]5[p-P«(^)] 

(18) 



The expression for the entropy 5„ has been given in Section III, Eqs. (9) and (13). 
Factors that do not depend on the overlap parameter Q^^ were omitted from Eq. (18). 

The hats on the functions Qa^ir,r') and p^(r) indicate that they are functions of the 
microscopic state {r^^it)} of the system. They are defined as in Eqs. (3) and (5), but 
without averaging: Q^p(r,r') = !dt8(r-r^(t))6(r'-rp(t)), and p^(r) = ! dt6(r -r^(t)) . 
Note that Qa^(r,r') appearing in Eq. (18) is just a dummy integration variable. It is 
convenient to write Eq. (18) in the following form: 



[Z"] = p!2„p(?,F') exp ^-I^^^dQ,^(r,7') exp[-|^] 

_ B B _| B 



(19) 
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The definition for E„ is clear from comparison with Eq. (18): 



— ^ = - — 2^J JFiJF2JF3JF4 2„p(Fi,F2)!2ap(''3''i)^('i'''3)^(''2'^4) (20) 

In the following we will refer to 5„ as the "entropy," and to as the "energy," 

although strictly speaking this is not correct. The great advantage of the independent 
interaction model as compared to the sequence model is the relatively simple expression 
for the energy. We will first present a highly simplified qualitative analysis that will 
reveal the existence of a critical value d^=6 for the dimension d . We will show that for 

d >d^ freezing can only take place at the microscopic length scale, implying the validity 
of the random energy assumption, while for d <d^ the freezing is of a more complicated 
nature. Following Ref. 15 we assume that (p has a characteristic scale R» a , and can, 
therefore, be written as 



X, 



n-l 



R R 



(21) 



Note that due to translation invariance the function (pj, which has unit scale and unit 

value, can be written as a function of (n - 1) arguments. The exponent of the scaling 

relation follows from the normalization ldx(p^(x) = N of (p (this normalization 

condition can be found by integrating the density p^^ir) given in Eq. (9) over the volume 

of the system). The idea is to find the scaling relations for the entropy and the energy, and 
to see whether the free energy has a stable extreme. The scaling relation for the entropy 

was derived in Ref. 15 to be 5 oc R~^ . The scaling relation for the energy Eq. (20) can be 

found by inserting for V(r) the generalized Coulomb potential V(r)o^ r^~^ , for the 

overlap parameter Q^^(x,y) the expression given in Eq. (9), and for (p(ii,---,x„) the 

scaling relation Eq. (21). After some work one finds eventually E oc r~'^'^ , and the free 
energy becomes 

For n<\ both constants Aj and are negative, and maximization of the free energy 
will yield a non-trivial value for i? if J < 6 (note that for short-range interactions the 
critical dimension is given by =2). In order to improve upon this rough analysis, it is 
necessary to use the more sophisticated trial function for cp that has been given in Eq. 
(10). First we rewrite expression Eq. (20) for the energy, using the translation invariance 
inside the polymer globule. 
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E„ Nf 1 



k^T 8p (iTif 



\dkV\k)Y^Ql^(k) (23) 



Inserting the trial function Eq. (12) for Q^^, and taking for V the screened Coulomb 
potential Eq. (17), one obtains 

_?^ = _«'-p^*^^^^^£ti^ (24) 

After taking the limit n — > the energy becomes 



k^T 



1 

= lim--^ = iVpK-V2VV8 {dxf{K^D(x)) 



/(0 = -V? +VV2 (1 + exp[V2] j^l-4>(V?/2) j (25) 

X 



4> is the error function. For future reference it is useful to give the series expansion and 
the asymptotic expansion^^ for / 

f(t) = -J%/2-2t^^^ for t«l 

(26) 

f{t) = 2r^l^ for t»\ 

Using the asymptotic expansion, an interesting observation can be made. It follows that 
for K D » 1 the energy can be approximated by 



1 





1 

^dxD'^l^ix) for Vd»k~^ (27) 



which coincides with the expression for the energy derived in Ref. 15 for the case of 
short-range interactions. Since according to Eq. (12) ^D^p is the freezing length scale 

(by which we mean the length scale at which the conformations of replicas a and [3 
coincide), one can conclude that with respect to the freezing behavior the interaction is 
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effectively short-range on length scales R»K~^ . Taking Eq. (13) for the entropy S 
and ignoring numerical factors of order unity, the final expression for the free energy F 
in the independent interaction approximation becomes 



-\^K(x) + c' \dxf{D(x)) 



NK^a^ •> X 



(28) 



c = Kp /k K{x) = x~^K(x) D(x) = x^D(x) 



Note that the dimensionless quantities K and D are still related via Eq. (14). If the 
independent interaction model is derived from the sequence model, then c is of order 
unity (Section V). However, we will analyze the model for general values of c , and find 
that non-trivial replica symmetry breaking occurs for c^l. The free energy F must be 
maximized, so a necessary condition for K is obtained by differentiating Eq. (28) with 
respect to K(x), using Eq. (14). By differentiating the result with respect to x, one 
obtains 



D(x)={f^''y 



^ K\x) ^ 
2c 



or 



dK 
dx 



= 



(29) 



where /'^^^ stands for the second derivative of / . Differentiating again with respect to x 
leads to the following solution to the maximization condition of Eq. (28): 



/ 



(3) 



2 

V J 



= -|K^x 

4 



or 



dx 



= 



(30) 



The reduced quantities x and K are defined by 



X = cx 



r\ = — = 

C C 



(31) 



The non-trivial solution given in Eq. (30) is shown graphically in Fig. 2. It exists only for 
X>Xg =1.49, and it has two branches. Since the function q(x) corresponding to the 

overlap matrix q^^ given in Eq. (12) has to be an increasing function of x , it follows that 

the upper branch is the physical one. The asymptotic behavior of K for large X can be 
found using the series expansion of / , leading to 



K = X for X ^ oo 



(32) 
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In terms of the original variable x, this solution exists on the interval X^/c<x<l, 
provided that c> x^. . However, one can check that the two expressions Eqs. (14) and (29) 
for D(x) do not coincide. Therefore, there rnust^^ be an interval Xq < x < 1 in which the 
second solution given in Eq. (30) is valid, that is, K is constant on this interval. The 
value of this constant Kq is, given X^/c <Xq <l, fixed by Eqs. (14) and (29). The 
maximization condition for the free energy, together with the physical condition that 
K(x) must be a non-decreasing function of x , imply that K(x) is continuous at x = x^, 

which fixes both jcq and Kq . It follows that = k"^^q must be a root of the following 
equation: 

/(2)(2/f )=F/2c^ (33) 

~3+i 

Eq. (26) implies that the left hand side of Eq. (33) is proportional to ^ for small K , 

while it is proportional to K^^^ for large K . This means that for small c Eq. (33) has no 
solution at all, while for large c it has two solutions, the larger of the two being the 
physical one. Using the series expansion for / , we obtain for c » 1 

K,=cy2 (34) 

Under the same condition c » 1 , the non-constant part of the solution follows from the 
combination of Eqs. (31) and (32) 

K(x) = c^x x>x^oc (35) 

The value for Xq follows from continuity at x = Xq: 

Xq=1/2 (36) 
The plateau value Dq for D is given by (use Eq. (14)) 

'''' 

Since x = l corresponds to self-overlap, it follows from Eqs. (12) and (37) that each 
replica is frozen on the self- screening length scale. The presence of a plateau means that 
there is a non-zero probability that two different replicas have perfect overlap down to the 
self-screening length scale. The non-constant part of K(x) , which is absent in the case of 

short-range interactions, corresponds to the possibility that the conformations of two 

replicas coincide down to length scale R , with K"^ <R< . The existence of this 
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interval for R is physically reasonable, because at length scales shorter than k" the 
thermal energy k^T is larger than the electrostatic energy, while at length scales larger 

than the interaction is effectively short-range due to the screening. The resulting 
shape of K(x) is shown in Fig. 3. Note that in the limiting case K = the overlap 
parameter q(x) increases continuously from zero, meaning that there are replica pairs 
having an arbitrarily small overlap. This implies that there is no upper bound for the 
"tube diameter" Z)~p (see Fig. 1), and one could describe the situation by saying that the 
freezing occurs on all length scales larger than the self-screening length. 

Let the overlap q^^ between two replicas a and P be defined as in Eq. (12). Let P(q) 
denote the probability density that two randomly picked replicas have overlap q . Then 
P(q) can be obtained^^ from the inverse of q(x) via P(q) = dx/dq. The dependence of 
P(q) on q is shown in Fig. 4. It has two delta peaks, one corresponding to the non-zero 

probability to find two replicas that have perfect overlap on the length scale k"^ , the 

other corresponding to the probability that two replicas have no overlap at all. The 
continuous part, which is due to the long-range character of the interaction, shows that 
there is a multitude of weakly overlapping states (compare with Ref. 15). 

V. SEQUENCE MODEL 

As explained in the previous section, self- screening does not exist in a system of particles 
interacting via independent coupling constants 5^ . Since at length scales R > k"' self- 
screening is very important in the polyampholyte globule, it is necessary to use a more 
realistic model in order to find the physical relevance of the results obtained in Section 
IV. In the sequence model the composition of the chain is described by its monomer 
sequence . The parameter i = l,---,N counts the charged monomers in the order 

determined by the chain connectivity. If 'd, =-i-l (1^, =-!) then the charged 
monomer has a positive (negative) electric charge of magnitude e . For computational 
convenience it will be assumed that i^, has a Gaussian distribution with mean zero and 

variance unity, i.e. the parameters -d, are drawn independently from the distribution 



P('&) = 




This simplification can be justified in the following way. On physical grounds, freezing 
might only occur at length scales R ^ K"^ , which is confirmed by the analysis in Section 

IV. Since a « K"' , it is possible to perform a renormalization transformation by defining 
coarse grained monomers each consisting of N-^ monomers, where 1 « « N^^ , 
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without changing the qualitative behavior of the model. The coarse grained monomers 
have a Gaussian charge distribution, which justifies the above-mentioned simplification. 
In the continuous representation of the chain^^ the partition function takes on the form 



(39) 



Again, it is assumed implicitly that the integration over chain conformations is restricted 
to conformations for which the density as defined in Eq. (5) is constant. For the potential 
V(r,r')=V(\r-f'\) we substitute the unscreened Coulomb potential V{r) = l/r. In 
order to calculate the free energy using replica theory^^ one has to average the partition 
function Z" of n non-interacting replicas over all possible sequences 1^(0 weighted 
according to the distribution given in Eq. (38). Following Refs. 16-18, it is convenient to 
introduce the parameter \jif(r) describing the charge density fluctuations: 



(40) 



It is defined as the deviation of the charge density from its average value (which is zero 
in our case). Since opposite charges attract each other (in Ref. 16 this situation was 
referred to as "the mixing case") the thermodynamic average of this parameter will 
always be zero. After writing the partition function in terms of \j/(r) , it is possible to 
make the effective Hamiltonian linear in ij by performing the Hubbard-Stratonovich 
trick, enabling the calculation of the average over the sequences. One obtains (the 
entropy 5„ has been defined in Eq. (18); see also Refs. 16,17 for additional details) 



[Z"] = |j^„(F)|rfv|/„(F) exp -^J^|jFJF'y(F,F')v|/„(F)v|/„(F') +/ |JF^„(F)V|/„(F': 



* f J(2„p(F,F')exp (drdr'Q,^(r,nK(r)%(n-^Y. Ur<^l(r) 

(41) 

This is a suitable moment to discuss the validity of the various approximations made in 
arriving at Eq. (41). One of the approximations is the assumption of a Gaussian 
probability distribution for the monomer charges, which causes the absence of terms of 

higher order in (|)„(F). If the partition function [Z"] is calculated using the discrete 
distribution ^(1^ = -l-l) = P(i^ = -1) = 1/2 , fourth and higher order terms are generated. As 
discussed earlier, the use of a Gaussian probability distribution can be justified by 
considering a coarse graining transformation. The higher order terms missed in this way 
are unimportant due to the fact that the thermodjTiamic average of ^a(^) is zero. 
Another approximation is the neglect of the fluctuations in the total density 
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p(r) = p+(r) + p_(r) , where p+(r) and p_(r) stand for the local density of positively, 
resp. negatively charged monomers. Note that it would not be easy to get rid of this 
approximation, since the trial function Eq. (10) does not incorporate density fluctuations. 
At first sight one may wonder whether it is consistent to take the fluctuations in the 
charge density \|/(r) = p_^(r)-p_(r) into account, while neglecting the fluctuations in 

the total density. However, the charge density and the total density play completely 
different roles. The fluctuations 5\|/ in the charge density are responsible for the 
electrostatic attraction^^ (note that 5\|/ = \|/); in other words the energy of a globule in the 
coarse grained state {\|/(r),p(r)} is a function of 5\|/(r) alone, and does not depend on 
5p(r) (see the first contribution to the effective Hamiltonian in Eq. (41)). Since this 

electrostatic interaction is responsible for the freezing transition, it is clear that the 
fluctuations 5\/(r) in the charge density are essential, and must be taken into account 

even in the lowest order approximation. The fluctuations 5p(r) in the total density, on 
the other hand, are far less important. As can be seen from Eq. (41), they appear only in 
the combination (p + 5p), for instance as argument of the entropy 'S'„((2otp>P(x)- At the 

relevant length scales R > k"^ , the fluctuations 5p are small compared to the average 
density p , and in the lowest order approximation it is justified to neglect them. Returning 
to Eq. (41), it is straightforward to proceed by calculating the Gaussian integrals over 
\|/p((r) and ([)«(?). The result is^^ 



[Z"]cx|rf2^p(F,F')exp 



H 



eff 



(42) 



^" ^ - Ji^^kidetA„p(^) 



k^T 2 (271)' 



The elements of the Parisi matrix A are given by 



AxpW = -eapW 



2 liYik) 



(43) 



Factors that are independent of the overlap parameter were omitted from Eq. (42). Since 
the chain is in a globular state, which is a well-defined state with short-range 
fluctuations, ' it is reasonable to assume that the saddle point method gives a good 
approximation to the integral over Q^^ (k) . Consequently, the equilibrium value of the 

overlap parameter can be obtained by finding a stable extreme of the effective 
Hamiltonian Hf^ . In the Umit n — > a stable extreme corresponds to a local maximum^^ 
of the mean- field free energy F given by 
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F = Mi-Hf (44) 

In order to calculate F , note that the expression for the entropy S is the same as it was 
for the independent interaction model; see Eq. (13). To find the expression for the energy, 
we use the following identit)^^ for Parisi matrices A 



lim llndetA.p(^) = ln(a-<a»+^- f^lnf l-^^l 
«->o n ^ a -{a) ix I a-\a) 



1 X 

~a=A^ <a> = |^a(x) [a](x) = xa(x)-|^ya(y) 



(45) 



The quantities a , {a) , and [a](x) all depend on A: . It is convenient to divide all matrix 
elements A^^p by p/2 ; from Eqs. (42) and (45) it follows that this corresponds to adding 

a constant to the free energy. Combining Eqs. (12), (43) and (45), and substituting for 
V{r) the unscreened Coulomb potential y(r) = 1/r , one arrives at 

a=l + ^- = l + -K~J-k^ a(x) = exp[-iD(x)A:^l 



Kp = 47ip^ 



(46) 



where Kp^ is the polymer self- screening length (note that p differs by a factor 2 from the 

charge density n used in Refs. 1,10). It was mentioned in Section IV that if the 
expression for the energy in the sequence model is expanded in powers of the potential, 
then the first term of this expansion gives the expression for the energy in the 
independent interaction model. Since Eq. (45) for the energy depends on the potential 
V{k) via the parameter a , an expansion in powers of the potential is equivalent to an 

expansion in powers of Ija . This expansion is justified for k -vectors with length 
k » Kp^ (note that since D{x) is positive, we have the inequalities < {a) < 1 , and 
0<a(jc)<l); corresponding to length scales shorter than the self-screening length. 
Under these conditions the energy Eq. (45) can be approximated by 

^^ = \aa+-^\dxa\x) if a»\ (47) 

kj 25^ J ' ' 

After neglecting terms that are independent of the overlap parameter, this expression can 
be rewritten as 
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k^T 2 (271)^ J 2(1 + x/k'f n^on^^ ' 

(48) 



Kp + A: ) is the Fourier transform of the screened Coulomb 

interaction with inverse screening length Kp. As mentioned, this expression for the 

energy coincides with the one derived from the independent interaction model (see Eq. 

(23)), with the important difference that the parameter does not stand for the salt 

screening length, but for the polymer self- screening length. Therefore, is not a free 
parameter: the globule density and the temperature fix its value. Following the same lines 
as in Section IV, one arrives at Eq. (28) for the free energy with the parameter c given by 

3 1 



[dkV\k-K)Yim-YQl^{k) 



If one inserts the expression for p as given in Eq. (1), one finds c = 0.5 . According to 
Section IV there exists a critical value =1.5 for c such that for c <c^ the non-trivial 
solution given in Eq. (30) does not exist (it can be shown that the critical value c'^ for the 
existence of a 1-step replica symmetry breaking solution has roughly the same value). Of 
course, the expression p-f^i/b"^ for the density is very crude: it is based on hand 
waving arguments and all numerical factors of order unity have been omitted. 
Nevertheless, due to the very weak dependence c <x p~^/^ of c on p we conclude that the 
present analysis predicts the absence of replica symmetry breaking in a weakly charged 
polyampholyte globule under -solvent conditions. Under poor-solvent conditions the 
situation is even worse since the density appears in the denominator of c . At first sight, 
this seems a bit contradictory: under poor solvent conditions the density is higher, and the 
electrostatic interactions are stronger. However, since the monomers are confined to a 
smaller volume, the constraints imposed by the polymer bonds become less important 
(the monomers are "asymptotically free"). They can rearrange themselves more easily in 
order to avoid unfavorable interactions, and the frustration is less than under - 
conditions. This reasoning can be made more quantitative by noting that under poor- 
solvent conditions the blob size is smaller than the self-screening length k"' : their 

ratio scales with the density p according to i?i,/i^p^ °^ • 

The expansion of E{k), leading to the independent interaction approximation to the 
energy, is justified only for k » k"' , which corresponds to length scales shorter than the 
polymer self-screening length. In order to study the freezing behavior on larger length 
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scales, one must resort to the full energy expression given by Eqs. (42), (43) and (45). By 
writing the entropy in terms of D(x) using 



K(x) 



X 

xD(x) + ^dyD(y) 



(50) 



and differentiating the free energy with respect to D(x) , it is possible to derive the 
following integral equation for D : 



X 

J 



dy 



D'(y) 



" bD(y) + ^dzD(z)V 



4% 



X oo 

dyD (y)\dpp (51) 



However, this equation is difficult to handle, even numerically, because for each x it 
involves the value of D(y) over the entire interval y& [0,1], making it impossible to 

grow the solution from x = . Therefore, instead of solving Eq. (51), we will show that it 
does not have a solution corresponding to an overlap parameter q(x) that increases 

continuously from zero. This means that there exists a maximum value for the "tube 

diameter" D^^ (see Fig. 1 for illustration), and there is no such phenomenon as "freezing 

on all length scales." To proof this point, first note that the second term on the right hand 
side of Eq. (45) is zero, because according to Eq. (13) we have K(0) = , so D(0) = oo , 

and a(0) = . By combining Eqs. (13), (42) and (45), one can write the expression for the 
mean-field free energy F = lim„^Q in the following way: 



4F 



1 



-^K(x) + - 



X 



3/2 



y'a' J 



dk k' 



(a) 
l + k 



H 



dx 



X 



[a](x) 



1 + A: -<a) J 



a(x) = exp[-^D(x)^^] 



k = K'^k 



K(x) = k;'K(x) D(x) = kID(x) 

(52) 

As a reminder we note that the first contribution to the free energy has an entropic origin, 
while the second contribution has an energetic origin. In order to arrive at this expression, 
the 2p(|3 -independent term ]na has been subtracted from E{k)lk^T , transforming the 

first term on the right hand side of Eq. (45) into ln(l-(a)/a ) . Note that D and K are 
still related to each other via Eq. (14). Substituting Eq. (1) for the density it follows that 
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the numerical factor in front of the energy is of order unity. As mentioned before, an 
important question is whether or not the sequence model, like the independent interaction 
model, permits a solution K(x) that increases continuously from zero. Physically this 
would mean that there are replica pairs having an arbitrary small, but non-zero overlap. 
The existence of such a solution is suggested by the fact that the unscreened Coulomb 

interaction has no intrinsic length scale. Since D = K^D and ^ = Kp ^ are rescaled 
quantities, overlap at a length scale which is large compared to the self-screening length 
corresponds to D » 1 , and so ^ « 1 . Since q(x) has to be an increasing function of 
X , the same must be true for K(x) . Also, according to Ref. 30 it is to be expected that 
there is an interval [xq,1] of x -values where K(x) is constant, and of order unity. 
Therefore, we look for a solution K(x) that increases continuously from zero at the left 
side of the interval [0, 1] , and reaches a plateau value ~ ^ the right side (the 
solution Fig. 3 to the independent interaction model has these properties if K = ). The 
most obvious thing to do is to assume that for small x -values the function K(x) scales 

like K(x) x""'^ . We must have m > 1 , and for m < 3 the dominant contribution to the 
integral representing the entropy comes from small x -values. Before inserting this type 
of solution into Eq. (52) for the free energy it is useful to make the following 
observations: 

• For the replica symmetric solution K(x) = one has a(x) = 0, and so both the 

entropy and the energy are zero. Therefore, in case of broken replica symmetry, the 
free energy must be positive. 

• The entropy, as given in Eq. (13), is always negative, and so is the first contribution 
to the energy. 

• The second contribution E2 to the energy in Eq. (52) is the only positive term, and 
must at least counterbalance the entropy in case of replica symmetry breaking. Note 
that E2 is smaller than the term £'2 obtained by omitting the -term from the 

denominator (Eq. (53)). Since this -term originates from the potential, removing 
this term corresponds to considering an infinitely strong interaction. 

Ei=M'UpJl-^^] S = -hKU) (53) 

J X J l-\a> J J X 

00^^ 

In the appendix it is shown that for functions K(x) which scale like x'"~^ for small x , 
the entropy S is dominant over E2 , which implies that it is dominant over the energy E , 
and the free energy as given in Eq. (52) is negative. Since the replica- symmetric solution 
has free energy zero, solutions which increase continuously from zero following a power 
law cannot exist. We believe that it is justified to make the stronger statement that the 
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maximization condition of the free energy given in Eq. (52) does not have a solution 
corresponding to an overlap parameter q(x) that increases continuously from zero. Note 

that since in the appendix it is proven that S is dominant over £'2 , our statement 
concerning the non-existence of such a solution applies to all situations for which 
V(k) >0 \/k. This condition is fulfilled if unlike particles attract each other ("mixing 

case"), and the (isotropic) interaction strength decays not slower than 1/r with the 

distance r . Although the exact solution to the maximization condition of the free energy 

may still be non-trivial, it follows from our analysis that the value of K(x) must be of 

order unity wherever it is non-zero. Therefore, it is not too bad an approximation to 
consider a 1-step solution, i.e. 

K(x) = x<x. 

(54) 

K(x) = Kq x> Xq 



Investigating this solution numerically shows that the factor in front of the energy (see 
Eq. (52)) should at least equal 8 in order to make the solution stable. It follows that in 
order to get replica symmetry breaking the density should be much lower than the one 
given in Eq. (1). If the chain is present in a poor solvent, the situation is even worse; see 
the discussion following Eq. (49). Therefore, we believe that it is safe to conclude that 
there is no replica symmetry breaking under these conditions. In the summary we discuss 
under what conditions replica symmetry breaking can occur, and draw some qualitative 
conclusions from the results obtained in the present section. 



VI. SUMMARY 

In this paper we investigated the possibility of freezing in the globular state of a 
"statistically neutral" random polyampholyte chain in a 6 -solvent for the situation where 
the polymer self- screening length is much larger than the average distance between two 
charged monomers (weakly charged case). Under these conditions the globule is far from 
being maximally compact and still contains a lot of solvent. Due to some general 
considerations (see the introduction) it was to be expected that the long-range character 
of the interaction changes the nature of the freezing transition. Using the analytical tools 
developed in Ref 15, we first studied the globule in the independent interaction 

approximation. This lead to the prediction that if the salt screening length is larger 

than the polymer self- screening length k"^ , the repUca symmetry breaking is non-trivial 

on length scales R satisfying k"^ <R< . However, one should be cautious in 

accepting the predictions made using the independent interaction model, because in this 
model self- screening does not exist. Therefore, we checked the results using the more 
realistic (and more complicated) sequence model. It was shown that on length scales 

smaller than K"^ the sequence model is equivalent to the independent interaction model, 
with the important footnote that the parameter , which represents the salt screening 
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length in the independent interaction model, now represents the self-screening length 
Kp^ . Adopting this physical interpretation of , the interesting regime K"^ <R< 

found in Section IV does not exist. In order to study the globule at length scales R^K.~^ 

one has to resort to the sequence model. We showed that for any interaction that favors 
mixing of unlike particles, and which decays not slower than 1/r , the overlap parameter 
q(x) is discontinuous in the sense that it must be of order unity wherever it is non-zero. 
Physically, this can be attributed to self- screening, which renders the interaction 
effectively short-range. For weakly charged polyampholytes in a -solvent or a good 
solvent we arrived at the conclusion that although the system is frustrated to a certain 
extent, no replica symmetry breaking occurs. In a poor solvent, the situation is even more 
relaxed: the frustration due to the presence of the polymer bonds is less if the monomers 
are confined to a smaller volume. At the same time, this suggests a way to increase the 
frustration in order to attain a state with replica symmetry breaking: by forcing the chain 
to remain in a state with smaller density than the one given by Eq. (1). In such an 
extended state the constraints imposed by the polymer bonds are more important, and the 
frustration is larger. This situation can be achieved in the following way. If the faction / 

of charged monomers increases, the density becomes larger. For a certain value of / , the 

density reaches its maximum value p^^^^ - f /b^ corresponding to the compact state. If 

/ increases beyond this value, the electrostatic attractions become stronger, but the chain 

is not able to respond by further collapse. So, it has to stay in the "extended" state 
corresponding to the maximum density. Of course, Eq. (1) breaks down before this point 
is reached. Also, in order to obtain an estimate for the fraction of charged monomers 
needed to attain the compact state, one should take into account the fact that in a dense 
globule the value of the Bjerrum length £ is much larger than its value £ ~ O.Vnm in pure 
water. Summarizing, the only possibility to attain replica symmetry breaking in 
polyampholytes is by considering maximally compact globules with a large fraction of 
charged monomers. However, in this case the screening length is of the same order of 
magnitude as the monomer size b , and the analytical analysis presented in this article 
breaks down. Nevertheless, we believe that the results of Section V remain qualitatively 
valid in the sense that the freezing transition can at least qualitatively be described by a 1- 
step replica symmetry breaking,^"^ and the system behaves in the same way as a globule 
with short-range interactions. Since we attribute this to the presence of self- screening, 
this conclusion is probably only true in the "mixing case", where unlike particles attract 
each other. In the segregating case, screening does not exist and the long-range character 
will have an important effect on the nature of the freezing transition. 
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APPENDIX A 



In this appendix we proof that if the function K(x) follows a power law K(x) <^ x'"~^ for 
small X , then the entropy S is dominant over the energy £'2 (the definitions for S and 
£'2 are reproduced in equation (Al)). As explained in the main text, the implication of 
this result is that in polymer globules with two kinds of monomers where unlike 
monomers attract each other, the overlap parameter q(x) must be of order unity 
wherever it is non-zero, even if the interaction is long-range (sufficient condition: 
V(k) > \/k). This means that there exists a maximum value for the "tube diameter" 
(see Fig. 1). 



dx 



X 



K(x) 



(Al) 



00 

=1' 







I -(a) 



We proceed as follows. For K(x) ^ we show that E(x) x'" ^ for small x , and 
find m' in terms of m . In order for the (positive) energy to be dominant over the 
(negative) entropy, m' must be smaller than m . We will show that this condition cannot 
be satisfied. Using the expressions for (a) and for [a](x) given in Eqs. (45) and (46), the 

expression for E(k) becomes (the tilde on the integration variable k will be omitted, and 

the prefactor c in D(x) = ex'"' is absorbed in A: ) 



00 

Eix) = ^dkk^\n 




^ xexp[-x-"k^]-n-^k^^T[-n-\x-''k^] ^ 
l-n-'k^f"r[-n-\k^] 



00 

r(s,x) = ^dy y'~^ exp[-3;] 



(A2) 



In order to find an upper bound for E(x), we construct a lower bound f(k) for the 
denominator in Eq. (A2); that is, we construct a function f(k) such that 



f(k)<g(k) = l-n-^k^^"r[-n-\k^] \/ k>0 (A3) 

The dependence of g(k) on k is shown in Fig. 5 for n = 2 , in which case it starts at zero 
with finite slope, and converges to unity for large k . For general values of n , the first 
term of the series expansion is 
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g(k) = -n-^r[-n-^]k^^'' + 



(A4) 



Note that r[-n ^ ] is negative. We split the interval 0<k<°° in two parts, separated by 

= x"^^. In the left part we choose f(k) to be given by the first term of the series 
expansion of g(k) , in the right part we choose f(k) to be constant. More precisely (see 
Fig. 5 for illustration for the case n = 2 for which f(k) is linear on < A: < )• 



f(k) = -rn~^T(-n-^)k^^" 
J{k) = -rn-'n-n-')x 



k < x"'^ 



(A5) 



r = 1-e 



e« 1 



Since the numerator in Eq. (A2) is a decreasing function of k for A; > x"^^ , the argument 
of the logarithm will not become negative. Substituting f(k) for g(k) in Eq. (A2) leads 

to the sum of two integrals. In both integrals, the integrand is a function of p = k x~"^^ , 
and we obtain 



\E(x)\<x^"'^^dpp^\n 



f -lln r 2t -It-t -1 2^^ 

p I exp[-/7 ]-w T[-n ,p ] 



+ 



+x 



3n/2 



OO 



p^ln 



exp[-/?^] - p^/"r[-n ~\ p^ ] 
- rn~^r[—n~^] 



-1 2-,\ 



(A6) 



For given value of n , the integrals in equation 1 are just numerical factors of order unity. 
Both terms in Eq. (A6) scale as x^"^^ . Since K(x)°'^x"~^, the negative entropy is 
dominant over the positive energy unless n<-2 , which is impossible since K{x) has to 
be an increasing function of x . 
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Legends to the Figures 



FIG. 1. Pictorial way to describe the trial function Eq. (12) for the overlap parameter. The 
conformations of replicas a and P coincide when the monomer positions are measured 

1/2 

with accuracy D^^ . When observed at better resolution, the conformations turn out to be 
different. 

FIG. 2. Qualitative representation of the non-trivial solution in Eq. (27). The rescaled 
quantities x and K are defined in Eq. (28). The solution has two branches. The dotted 
branch must be rejected since it leads to an overlap parameter q(x) that decreases with 

increasing x. The solid branch represents the physical solution; it converges to the 
asymptot K = x (dashed Une). 

FIG. 3. Approximate shape of K{x) calculated within the independent interaction 
approximation for the situation k~' » k"' . The jump in K is due to the absence of 

freezing at length scales larger than K~' . The interval where K is linear corresponds to 
length scales in between K~' and K~' ; in this range the freezing is qualitatively different 
from that occurring in polymer globules with short-range interactions. The plateau on the 
right corresponds to freezing on the length scale . If the independent interaction 

model is obtained from the sequence model by expansion of the free energy in powers of 
the potential, then k = Kp , and the solution shown here does not exist for weakly charged 

polyamphol5^es. 

FIG. 4. This figure shows, for the independent interaction model, the probability density 
P that the overlap between two randomly picked replicas equals q . The delta peak at the 
left represents the non-zero probability that two replicas have no overlap at all. The delta 
peak at the right represents the probability that two replicas have perfect overlap down to 
the self- screening length scale. The probability to find a larger overlap is zero, because at 
length scales smaller than the self- screening length the thermal energy prevents freezing. 
The continuous part of P(q) represents non-trivial freezing due to the long-range 
character of the interaction (compare with Ref. 15). 

FIG. 5. The solid curve represents the denominator g(k) = (I - n~^k^^"r[-n~^ ,k^]) of the 
quotient in Eq. (A2) for the case n = 2; the dashed curve represents its lower bound 
f(k) used to find an upper bound for the energy E(x) . For n^2 the dashed curve will 

not be linear for small k but increase according to the power law f(k)<xk^^". 
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